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INFRARED CLASSIFICATION OF GALACTIC OBJECTS 

Zeljko Ivezic^ and Moshe Elitzur^ 

(3 ABSTRACT 
O _ 

^ . Unbiased analysis shows that IRAS data reliably differentiate between 

5_i ! the early and late stages of stellar evolution because objects at these stages 

^ clearly segregate in infrared color-color diagrams. Structure in these diagrams 

is primarily controlled by the density distribution of circumstellar dust. The 
^ ! density profile around older objects is the steepest, declining as r~^, while 

young objects have profiles that vary as r~^/^ and flatter. The different density 



> 



profiles reflect the different dynamics that govern the different environments. 



, Our analysis also shows that high mass star formation is strongly concentrated 

OO 

^ ! within ~ 5 kpc around the Galactic center, in support of other studies. 

m 
o . 

O . Subject headings: catalogs — stars: AGB and post-AGB — stars: formation 



infrared: general — Galaxy: stellar content — Galaxy: structure 



O 

2 ■ 1. INTRODUCTION 

^ : The IRAS all-sky survey provides a unique opportunity to classify the infrared 

properties of astronomical objects from a homogeneous data set obtained with a single 
^ '. facility. Cross-correlations of various catalogues with the IRAS Point Source Catalogue 

(PSC) showed that certain Galactic objects tend to cluster in well defined regions of IRAS 
color-color diagrams. Notable examples include HII regions (e.g. Hughes & MacLeod 1989; 
Wood & Churchwell 1989 [WC]) and AGB stars (van der Veen & Habing 1988; VH). The 
reason for this clustering was not understood, nor was it clear whether such biased analysis 
based on pre-selection implies reliable selection criteria. From detailed modeling of dusty 
winds we were able to validate the VH selection criterion proposed for AGB stars and 
to explain its origin (Ivezic & Elitzur 1995; IE95). Here we extend this approach to all 
Galactic PSC sources in an unbiased analysis of IRAS fluxes without prior selections. 
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2. IRAS DATA AND ITS CLASSES 

There are 6338 PSC sources with flux quahty of at least 2 in all four IRAS channels. 
Removing the 651 sources identified as extra-galactic (Beichman et al 1985) produces our 
basic data set of 5687 Galactic IR objects. For AGB stars we have shown (IE95) that 
a high IRAS quality docs not guarantee the flux is intrinsic to the point source itself; 
high-quality 60 and 100 fim PSC fluxes can have a cirrus origin instead. Figure 1 shows 
that the problem afflicts all sources, not just AGB stars. Its x-axis is cirrS/Fgo, a measure 
of the ratio of cirrus background noise to the 60 fim signal. Intrinsic fluxes should have 
nothing to do with background emission, yet the [100-60] color is strongly correlated with 
cirr3/F6o when this noise indicator exceeds a certain threshold. The PSC fluxes of sources 
above this threshold reflect cirrus, not intrinsic emission. Following IE95, we remove all 
sources with cirr3 > 2 Fqq to eliminate cirrus contamination. This leaves 1493 objects that 
can be considered a reliable representation of Galactic infrared point sources. 

We submitted these data to the program AutoClass for an unbiased search of possible 
structure. AutoClass employs Bayesian probability analysis to automatically abstract a 
given data base into classes (Goebel et al 1989). AutoClass detected that the PSC sources 
belong to four distinct classes that occupy separate regions in the 4-dimensional space 
spanned by IRAS fluxes. Figure 2a shows the 2-dimensional projection of this space onto 
the [100-60]-[25-12] color plane. The IRAS colors of black body emission, which are 
independent of the black body temperature as long as it exceeds ~ 700 K, are marked by 
the large dot. There is hardly any data at that point. Instead, the data are spread far away 
from it, indicating that IRAS fluxes are dominated by surrounding dust. But rather than 
random scatter, the data show clear structure, with the four AutoClass classes occupying 
well defined color regions. As a further check on the reality of these classes we constructed 
their Galactic distributions, shown in Figure 2b. 

Classes A and B clearly separate in the color-color diagram. Classes C and D, on 
the other hand, are mixed together and are distinguished by their flux levels — class C 
fluxes are typically 10-100 times higher in all 4 bands. In principle, a single family could 
produce such behavior if split into nearby and distant objects, but this is not the case here. 
Extragalactic and heliocentric selection effects arc ruled out by the Galactic distributions. 
Both classes are comprised of Galactic disk objects, but the median flux of class C is 16 
times higher. If the two were drawn from the same population, class C sources would be on 
the average 4 times closer and their Galactic latitude distribution 4 times wider. Instead, 
the latitude histograms are essentially the same, both are centered on 6 = 0° with full-width 
at half-maximum 1.5° ± 0.3° for class C and 1.5° ± 0.1° for class D. Class D fluxes vary 
by more than 10^, consistent with a population distributed throughout the entire Galactic 
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plane.0 In contrast, class C fluxes vary by a factor of only ~ 20 and exceed 10^ Jy at 
both 60 and 100 /im for all sources. This is consistent only with a population confined 
to, at most, ~ 5 kpc around the Galactic center with a minimal luminosity of 10^ Lq. 
The Galactic longitude distributions corroborate these conclusions. Although both class C 
and D sources channel most of their luminosity to the far-IR, their spatial and luminosity 
distributions are very different. 

3. THEORY 

In IE95 we identified the cause of the particular IRAS colors of AGB stars: The dusty 
wind problem possesses general scaling properties. For a given dust composition, the 
solutions are predominantly determined by a single input parameter — the overall optical 
depth. All other input is largely irrelevant. As a result, the solutions occupy well defined 
regions in color-color diagrams. Indeed, the IRAS measured colors of AGB stars fall in the 
regions outlined by wind solutions with visual optical depth Ty ^ 100. 

In Ivezic & Elitzur (1997; IE97) we showed that scaling is a general property of 
radiatively heated dust under all circumstances, not just in AGB winds. The most general 
dust radiative transfer problem contains only two input quantities whose magnitudes 
matter — Ty and the dust temperature Ti at some point. All other input is defined by 
dimensionless, normalized profiles that describe (1) the spectral shape of the external 
radiation, (2) the spectral shape of the dust absorption and scattering coefficients, and (3) 
the dust spatial distribution. Physical dimensions such as, e.g., luminosity and linear sizes 
are irrelevant. Scaling applies to arbitrary geometries, and in IE97 we conducted extensive 
numerical studies of its consequences for spherical shells. With a black body spectral shape 
for the heating radiation, the temperature T hardly affects IRAS colors when T > 2000 K. 
The reason is that T is much higher than the Planck equivalent of the shortest IRAS 
wavelength (12 /im), which is only 240 K. Although this result was obtained in spherical 
solutions, its validity is general. Since neither its luminosity nor spectral shape matter, the 
heating source is quite irrelevant for the IRAS colors. Different objects segregate in IRAS 
color-color diagrams not because they have different central sources, but rather because 
their dust shells are different. 

In our extensive study of spherical shells we examined also the effect of Ti, which was 
selected as the dust temperature on the shell inner boundary. As with the heating radiation, 
and for the same reasons, Ti hardly affects IRAS colors when varied over the plausible range 
of dust sublimation temperatures 700-2000 K (IE97). This leaves the dust's composition, 
optical depth and spatial distribution as the only potentially significant reasons for the 
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variation of IRAS colors among different families of Galactic objects. However, dust 
properties generally do not show large variations among different objects and thus cannot 
be expected to induce substantial variations in infrared emission. The same applies to Ty, 
whose range of values is similar in most families. We conclude that objects segregate in 
IRAS color-color diagrams primarily because their spatial dust distributions are different. 

We corroborate these conclusions in detailed modeling with the code DUSTY (Ivezic, 
Nenkova & Ehtzur 1999). In these calculations, a point source emitting as a black body 
with temperature 5000 K is surrounded by dust whose absorption and scattering coefficients 
are those of standard interstellar mix. The dust is distributed in a spherical shell whose 
temperature on the inner boundary is Ti = 1000 K. The shell density profile is taken as 
a power law r~^, with p a free parameter. Scaling ensures that the solution set for each 
p is a one-parameter family characterized by Ty, leading to the distinct tracks shown in 
the [60-25]-[25-12] color-color diagram in figure 3. These tracks properly delineate the 
distribution of IRAS sources and we have verified that the same applies also to the [100-60] 
colors shown in fig. 2. Tracks for different p differ in the color-color regions that they cross 
and in the distance induced by Ty variation; the steeper the density profile, the smaller 
is the distance along the track for the same change in Ty. Both properties arise from the 
differences in relative amounts of material placed at different temperatures, and both are 
reflected in the data. Class A sources are well explained by the p = 2 track, indeed they 
cluster close to the black-body point. In contrast, classes C and D are explained by tracks 
of much flatter density distributions and are located rather far from the black-body colors; 
it takes only Ty ~ 0.1 to move an object along these tracks all the way from the black body 
point to the region populated by IRAS sources. 

Class B cannot be reasonably explained by the same models, even though its colors can 
be produced by extending the p — 2 track beyond Ty — 100. In addition to the implausible 
optical depths this requires, class B is unhkely to be the high-ry end of class A given 
their different Galactic distributions. An alternative to ry-variation is to modify colors by 
removal of hot dust. The displayed tracks have Ti = 1000 K, as appropriate for shells whose 
inner boundaries are controlled by dust sublimation, but in detached shells Ti is both lower 
and arbitrary. The dashed-line tracks in figure 3 show the effect of lowering Ti on the p 
= 2 track. The difference between Ti = 1000 K and 300 K is marginal because both are 
higher than the Planck equivalent of 12 /xm. However, further reduction in Ti alters the 
track significantly, adequately explaining class B colors. 

4. DISCUSSION 

Our modeling confirms that the primary reason for different IRAS classes is different 
dust density distributions. Only 5% of the sources require an additional variation of the 
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dust temperature on the shell inner boundary. Spherical shells were employed here as the 
simplest method to model extended, three dimensional dust distributions. Power laws were 
used just for illustration purposes and although their actual values should not be taken 
literally, they provide a good indication of the overall behavior of the density distributions 
in the four IRAS classes. The presence of disks, expected in various sources, should not 
significantly affect our conclusions. The standard disk spectrum F^, oc ly^^^ produces a single 
point in color-color diagrams and thus cannot effect the observed scatter. This spectrum 
is modified if the disk is embedded in an extended shell, but then the disk is expected 
to dominate only at sub-mm and mm wavelengths, longer than those observed by IRAS 
(Miroshnichenko et al 1999). Flared disk emission can be shown equivalent to that from a 
fiat disk embedded in an appropriate spherical shell. Therefore, such a configuration, too, 
cannot modify our conclusions. 

We have queried the SIMBAD database about our sample sources. SIMBAD 
identifications are occasionally ambiguous ("maser"), sometimes less than informative 
("part of cloud"), and reliability is not always certain. Nevertheless, they provide useful 
clues when there are major trends in the data. Class A returned the most decisive results 
— 88% of its members have possible optical identifications, of which roughly 90% are 
commensurate with AGB stars. Since class A obeys the VH criterion, this corroborates our 
earlier finding (IE95) that this criterion is both sufficient and necessary for AGB selection. 
In IE95 we present a thorough analysis of the IRAS fiuxes of AGB stars, including color 
tracks geared specifically for these objects (single chemistry grains, appropriate stellar 
temperature, etc.). While that analysis was necessary to verify the VH criterion, the p = 
2 track presented here captures the essence of the more detailed study, demonstrating that 
the density distribution is the leading factor in controlhng IRAS colors. These circumstellar 
shells have the steepest density distribution, setting them apart in the color-color diagrams. 
In addition, the color tracks are primarily controlled by a single parameter, Ty, hence the 
compactness of this class color region. 

Class B had the same identification rate but its composition not as homogeneous. 
Planetary nebulae comprise 40% of positive and possible identifications. At 13%, the only 
other significant group is "emission-line stars", a classification consistent with planetary 
and refiection nebulae. The remaining identifications span a variety of objects, indicating 
that detached shells may occasionally form under different circumstances. Significantly, the 
optical depths required for class B are lower than for the others, the dashed-line tracks in 
fig. 3 terminate at Ty = 10. 

SIMBAD identification rates for the two other classes are much lower — only 32% 
for class C and 21% for D. Among the 38 class C identifications, HII regions comprise 
the only significant group with 15 (40%). In class D, too, HII regions comprise the single 
largest group with 23% of the identifications, followed by young stellar objects at 11% 
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and planetary nebulae at 10%. These two classes are clearly dominated by star formation 
and early stages of stellar evolution, in agreement with their Galactic distributions and 
with previous attempts to associate star-forming regions with IRAS colors. In the most 
extensive study of this kind, 83% of 1302 IRAS selected sources were found to be embedded 
in molecular clouds and thus trace star formation, and the selection criterion ([25-12] > 0, 
[60-25] > 0.4) essentially identifies classes C and D (Wouterloot & Brand 1989; Wouterloot 
et al 1990). 

The IRAS colors of classes C and D imply dust density distributions flatter than for 
AGB stars. These colors are spread over large regions, reflecting variation in density profiles 
in addition to optical depth. The spread in [100-60] colors is smaller than in [60-25] because 
all shells are optically thin at both 100 and 60 //m while their 25 /xm optical depth can 
become significant (IE97). Class C colors occupy a sub-region of class D and are produced 
by the optically-thick end (ry ^ 1) of fiat density distributions. Objects whose colors fall 
in that region can belong to either class C or D. However, among all sources with class D 
colors, those with high fiuxes (> 1000 Jy at both 60 and 100 /xm) concentrate in a compact 
color region, hence the separate class C. Since all class C sources have L > 10^ Lq, they are 
high-mass objects and their concentration in the inner ~ 5 kpc of the Galaxy shows that 
the high-mass star formation rate decreases with distance from the Galactic center. This 
result is in agreement with studies of the initial mass function inside and outside the solar 
circle (Garmany et al 1982; Wouterloot et al 1995; Casassus et al 1999). 

Our results explain the findings of all earlier studies that were based on object 
pre-selection, and reveal the limitations of that approach. By example, consider the 
WC study. After identifying IRAS counterparts of known ultracompact HII regions, WC 
proposed the corresponding colors as a necessary and sufficient selection for all ultracompact 
HII regions and proceeded to estimate the birthrate of O stars. Codella et al (1994) then 
found that most HII regions in a more extended catalog indeed obeyed the WC color 
criterion. However, that included also diffuse, not just compact HII regions, therefore WC 
overestimated the O star birthrate. This clearly demonstrates the shortcomings of any 
classification based on a pre-selected population. The unbiased analysis presented here 
shows that IRAS colors refiect primarily the dust density profiles of circumstellar shells and 
provide a unique indication of the underlying object only in the case of AGB stars. IRAS 
data in itself is sufficient for differentiating young and old stellar objects; apart from a 
limited number of detached shells, IRAS sources belong to two distinct groups as is evident 
from both the color-color diagrams and the Galactic distributions: (1) class A sources are 
at the late stages of stellar evolution and (2) class C and D sources are objects at the early 
evolutionary stages. This differentiation occurs because the density distributions of dust 
around young stellar objects have flatter profiles, reflecting the different dynamics that 
govern the different environments. 
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Fig. 1. — [100-60] color vs. cirrS/Fgo, where Fqq is the 60 //m IRAS flux and cirr3 is the 
IRAS surface brightness at 100 //m around the point source. The flux quahty for all sources 
is at least 2 in all four IRAS bands, and 3 for those marked with dark dots. In spite of the 
high flux quality, the strong correlation with the cirrus emission at large values of cirrS/Feo 
is evident. 

Fig. 2. — (a) [100-60]-[25-12] color-color diagram for all the Galactic IRAS sources free from 
cirrus contamination. The four classes identified by AutoClass are marked with different 
symbols: Class A (282 sources) — crosses, B (80 sources) — circles, C (119 sources) — 
squares, D (1013 sources) — dots. The large dot marked BB denotes the colors of black 
body emission with temperature T > 700 K. (b) Galactic distributions of the four classes. 

Fig. 3. — [60-25]-[25-12] color-color diagram. Data are for the same sources presented in 
figure 2 using the same symbols. The full fines are model tracks for spherical dust shells 
with Ti — 1000 K on their inner boundary and density profiles with p — 0, 1, ^2 
and 2, as marked. Position along each track increases with dust optical depth ry away from 
the common origin ry = at the black-body colors, with ry = 0.1 marked by solid triangle 
and Ty = 10 by solid square. The cndpoint for each track is ry = 100. Dashcd-linc tracks 
show the effect on the p = 2 track of lowering Ti to 300 K (lowermost track), 200 K, 150 K 
and 100 K (uppermost track). These tracks end at ry = 10. 
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